Percolation on correlated networks 
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We reconsider the problem of percolation on an equilibrium random network with degree-degree 
correlations between nearest-neighboring vertices focusing on critical singularities at a percolation 
threshold. We obtain criteria for degree-degree correlations to be irrelevant for critical singularities. 
We present examples of networks in which assortative and disassortative mixing leads to unusual 
percolation properties and new critical exponents. 
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I. INTRODUCTION 

Real- world networks are correlated [l|, 0, H, 0, @ . Cor- 
relations between degrees of vertices in a network essen- 
tially characterise its structure. Various real-world net- 
works are markedly different in respect of degree-degree 
correlations [EIESBQ- In particular, social networks 
show assortative mixing, i.e., a preference of high-degree 
vertices to be connected to other high-degree vertices, 
while technological and biological networks are mostly 
disassortative, i.e., their high-degree vertices tend to be 
connected to low-degree ones [91- However, even the sim- 
plest correlated networks with only pair correlations be- 
tween the nearest-neighbor degrees are still poorly under- 
stood. Our aim is to find when the critical singularities 
for correlated networks of this kind coincide with those 
for well-studied uncorrelated networks and when and how 
much they differ. 

At the present time it is well established that the small- 
world effect and heterogeneity influence the cooperative 
dynamics and critical phenomena of models defined on 
the top of complex networks [l^. However, numerous 
studies were devoted mostly to more simple, uncorrelated 
networks. In an uncorrelated network with a heavy-tailed 
degree distribution, critical singularities of a continuous 
phase transition arc characterized by model dependent 
critical exponents which differ from the standard mean- 
field ones and the critical exponents of two and three 
dimensional lattices, see, for example, Refs. pdl [T3. fisj. 
The critical behavior depends on an asymptotic behav- 
ior of a degree distribution at large degrees. For per- 
colation on an uncorrelated complex network that was 
demonstrated in [13, [3 ■ One should expect, however, 
that for dynamical processes taking place in a complex 
networks, correlations are important. The simplest par- 
ticular kind of correlations in a networks are correlations 
between degrees of two nearest neighbors in a network — 
so called degree-degree correlations. In this work we con- 
sider only these specific, though representative, correla- 
tions. Inv estig ations of percolation 0, [l^ and epidemic 
spreading (TtT fisj demonstrated that the degree-degree 
correlations strongly influence these phenomena. The 
birth and growth of the giant connected component sig- 
nificantly depends on the type of correlations — whether 
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FIG. 1: Schematic representation of eigenvalues (crosses) of 
the branching matrix, Eq. (j4]|, in the limit N oo. There 
are four situations: (a) the largest eigenvalue Ai is nonzero 
while all other eigenvalues are zero, Ai = for i > 2; (b) both 
the largest eigenvalue Ai and the second largest eigenvalue A2 
are finite; (c) Ai diverges but A2 is finite; (d) the sequence 
of eigenvalues is unbounded. Case (a) takes place for uncor- 
related networks but also can occur in specific disassortative 
networks. Cases (b)-(d) occur only in correlated networks. 



the degree-degree correlations are assortative or disas- 
sortative. Compared to an uncorrelated network with 
the same degree distribution, the assortative correlations 
increase the resilience of a network against random dam- 
age, while the disassortative correlations diminish this 
resilience. 

One can construct an equilibrium network, where only 
degree-degree correlations are present — the maximally 
random network with given degree-degree correlations. 
This is impossible for non-equilibrium, in particular, 
growing networks. Growing networks necessarily demon- 
strate a wide spectrum of correlations and not only pair 
correlations between degrees of the nearest neighbors. 
This type of heterogeneity may result in an anomalous 
critical effect at the birth point of a giant connected 
component [H, [13, EH, [H, [2I, [H [H [21] . This transi- 
tion resembles the Berezinskii-Kosterlitz-Thouless phase 
transition in condensed matter [13, [ll] . The transition 
is explained by a specific large-scale inhomogeneity of 
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TABLE I: Critical behavior of the order parameter j/ oc (p — Pc) , where p is an occupation probability, in correlated networks 
with a degree distribution P{q) ~ Ai and A2 are the largest and second largest eigenvalues of the branching matrix, 

Eq. (j4]|. "Weakly" correlated networks, both assortative and disassortative, have the same critical behavior as uncorrelated 
network with the identical degree distribution, e.i., correlations are irrelevant. Model (144 |l with assortative mixing has a new 
critical exponent /3 in the range — 1 < a < 0. Here, the parameter a determines the average degree of the nearest neighbor of a 
vertex with q connections: q^-^^iq) = const + q~" ■ At a > this model is "weakly" correlated. Model (|62p with disassortative 
mixing has a new critical behavior at a > 0. For this model, q^^{q) = 1 + The (7, a) phase diagrams of models (I44|l and 
(|62p are shown in Figs. [5] and [H respectively. 
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non-equilibrium networks and related correlations. The 
large-scale inhomogeneity here means the difference in 
properties of vertices according to age [H, HO] ■ 

One should note that degree-degree correlations may 
also arise in an equilibrium network if self-loops and mul- 
tiple connections are forbidden. In particular, this was 
demonstrated in Ref. |3l||for the static model of scale- 
free networks [3^ . Noh[33j observed an unusual critical 
behavior in an exponential random graph model with a 
tunable degree-degree correlations. 

The earlier investigations B m, 113, m mostly focused 
on the effect of degree-degree correlations on the percola- 
tion and epidemic thresholds. In the present paper we in- 
vestigate the effect of degree-degree correlations in equi- 
librium networks on critical singularities at the percola- 
tion threshold. We demonstrate that the critical behav- 
ior is determined by the spectrum of a so-called branch- 
ing matrix and a degree distribution. The eigenvalues 
of this matrix are real and can be ordered in descending 
order. The largest eigenvalue determines the percolation 
threshold and can be both finite and infinite depending 
on a degree distribution in agreement with Refs. [ol. [l6j. 
We derive necessary and sufficient conditions for degree- 
degree correlations to be irrelevant for critical singular- 
ities. We give examples of strongly correlated networks 
with assortative and disassortative mixing in which at 
least one of these conditions is not fulfilled. These net- 
works demonstrate new critical singularities, see Table ID 
In particular, we propose analytically treatable assorta- 
tive networks with an unbounded sequence of eigenvalues 
in the infinite size limit and show that this peculiarity 
of the spectrum leads to a new critical behavior. Re- 
markably, this network may be robust against a random 
damage even if the second moment of its degree distribu- 
tion is finite. This is in contrast to uncorrelated networks 



which are robust only if the second moment of a degree 
distribution diverges. We also study specific disassorta- 
tive networks which demonstrate new critical singulari- 
ties and are fragile. Remarkably, this takes place despite 
the second moment of their degree distribution may be 
divergent and the largest eigenvalue is finite. Our results 
are summed up in Table [T] and Figs. and [5] which dis- 
play phase diagrams of two model networks with strong 
assortative and disassortative mixing. 

The paper is organized as follows. In Sec. |ll] we give 
a general description of correlations in a network and in- 
troduce the branching matrix. In Sec. Illll we introduce a 
basic equation describing percolation on a degree-degree 
correlated network and reconsider the effect of assorta- 
tive and disassortative mixing on the percolation thresh- 
old. In Sec. IIVI we find criteria for degree-degree corre- 
lations to be irrelevant for critical singularities. Sec. |V] 
introduces a simple model of degree-degree correlated 
networks which have a modular structure and an un- 
bounded sequence of eigenvalues of the branching ma- 
trix. In Sec. IVII we demonstrate that these networks 
have new critical singularities. Disassortative networks 
with unusual critical properties are studied in Sec. IVIII 
A detailed analysis of the spectrum, its relationships with 
structural coefficients of a correlated network and calcu- 
lation of the entropy of the model network are presented 
in Appendices El El and |Cl respectively. 

II. DEGREE-DEGREE CORRELATIONS 

We consider a random locally tree-like network of N 
vertices in which only pair correlations between nearest 
neighbor degrees q and q' are present. We assume that 
there are vertices with degrees q = 91,(727 ■■■qcut, and in 
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total there are Nd different values of degrees. This corre- 
lated network is completely described by a given symmet- 
ric joint degree-degree distribution, P{q,q') = P{q',q), 
otherwise the network is homogeneously random. The 
degree distribution P{q) can be calculated from P{q, q') 
as follows: 



Y^P{q\q)=qP{q)/{q), 



(1) 



It only depends on the degree distribution P{q). An in- 
tegral characteristic of degree-degree correlations is given 
by the Pearson coefficient : 



Y}B{q)-B]q{q-l)P{q), 



(9) 



where 



where the brackets (...) denote an average over the degree 
distribution P(g), for example, {q) — J2q1-Pil)- Below 
for brevity we use the notations: zi = (q) and Z2 = 
{q{q — 1)). It is convenient to introduce the conditional 
probability. 



ZlP{q,q' 
qP{q) 



(2) 



that if an end vertex of an edge has degree q, then the 
second end has degree q' . The functions P{q,q') and 
P{q'\q) are normalized: 



Y,P{q^q') = Y^P{q'\q)^l. 



(3) 



q,q' 



q' 



In uncorrelated networks the conditional probability does 
not depend on q: P{q'\q) — q' P{q') / zi. We define a non- 
symmetric branching matrix as follows: 



q q- 



(4) 



(5) 



Bqq. = [q' - l)P{q'\q). 

This matrix has the following property: 

q{q - l)P{q)Bqq, = q'{q' - l)P{q')B 

According to definition ([¥]), an entry Bqqi of this matrix 
is equal to the branching coefficient g' — 1 of an edge, 
which emanates from a vertex with degree q and has a 
vertex with degree q' at the second end, multiplied by the 
probability P{q'\q) that the second end has degree q' . In 
Sees, mil and IIVI we will show that the structure of the 
spectrum of this matrix determines the critical properties 
of the percolation. Relationship of this spectrum with 
clustering coefficients and the mean intervertex distance 
in a correlated network is considered in Appendix |Bl 

Using the branching matrix we can calculate the aver- 
age branching coefficient B{q) of an edge which emanates 
from a vertex of degree q: 



B{q) 



^Bqq'- 

q' 



(6) 



The coefficient B{q) is related to the average nearest- 
neighbor's degree of the vertices of degree q: 



(In 



M = Y.^'P{q\q) = B{q) 



1. 



A mean branching coefficient B is equal to 



(7) 



(8) 



qP{q) 



Zl 



qPjq) 

Zl 



(q) (q)' 



is for normalization. The Pierson coefficient is posi- 
tive for assortative mixing and negative for disassortative 
mixing. 

In an uncorrelated complex network, ^„„(<z) does not 
depend on q: qnniq) — B + 1. In contrast, in networks 
with assortative mixing, 5„„ (q) increases with increasing 
q while in disassortative networks it decreases. For exam- 
ple, qnniq) shows a power law decay qnniq) °^ q~'^'^ for 
the Internet Recursive scale- free networks, growing 

by the linear preferential attachment mechanism, have 
growing asymptotics 9„„((j') ~ Ing at large q at 7 > 3 
and decaying g„„ {q) as a power law at 2 < 7 < 3 [2l|, . 



III. TREE ANSATZ EQUATIONS AND 
PERCOLATION THRESHOLD 



A key quantity in the percolation problem is the prob- 
ability Xq that if an edge is attached to a vertex of degree 
q, then, following this edge to its second end, we will not 
appear in a giant connected component . The number 
of unknown order parameter components Xq is equal to 
the number of different degrees, Nd- In contrast, only a 
one-component order parameter x describes percolation 
on an uncorrelated complex network [sj |35| . 

Let p be the probability that a vertex is retained in a 
randomly damaged network. Within a tree-ansatz the- 
ory, which assumes that a network has a locally tree-like 
structure, equations for the probabilities Xq and the rel- 
ative size S" of a giant connected component have the 
following form [l6j : 



Xq = i-p + pJ2Piq'\q)i->^ 

q' 

i-s^i-p+pJ2Piq)i^q)' 



(10) 



(11) 



Equations (fTO|) and (|lip directly generalize equations de- 
rived for percolation on an uncorrelated complex net- 
work [13,1^, where Xq = x, see also Refs. 2, 3, 10]. 
The set of Eqs. (|10p determines unknown probabilities 



for q — qi, q2, ■■■qcut- Newman ^ originally derived 
these equations using generating functions, and numeri- 
cally solved them for various networks. The analysis of 
Eqs. (|10p shows that the birth of the giant connected 
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component is a continuous phase transition. Below the 
percolation threshold, i.e., at p < p^, these equations 
only have a trivial solution Xq — \, and there is no gi- 
ant connected component. A giant connected component 
is present above the percolation point, p > Pc, where 
there is a solution with < 1. Introducing a parameter 
j/g 1 — Xg, we rewrite Eq. (jlOp as follows: 

q' q'=3 n=2 ^ ^ 

(12) 

One can solve this set of equations near the critical point, 
T = p/pc — 1 -C 1 when < 1. 

First we study the percolation threshold for a degree- 
degree correlated network. We will use the spectral prop- 
erties of the branching matrix, Eq. (g]), so let us remind 

some basics. Eigenvalues Xi and eigenvectors $g ^ associ- 
ated with these eigenvalues are defined by the equation: 

J2Bqq'<i>':;}^\<^>';^\ (13) 

where the index i labels the eigenvalues, i — 1,2, . . . , Nd- 
We consider the case of positive entries: Bqgi > 0. The 
following statements hold, (i) The largest eigenvalue Ai is 
positive, (ii) The entries 'Pq^-' of the maximal eigenvector 

are positive, > 0. (iii) All eigenvalues Ai are real 
and can be ordered: Ai > A2 > ... > Xn^- (iv) The 
eigenvectors associated with these eigenvalues form 
a complete orthonormal basis set. That is, 

^«;((Z)$«$W=5,,„ (14) 

W{q)j2<f[^^^[^^ ^6q^q,, (15) 

A 

where w{q) — q{q — l)P{q)/zi is a weight function, 6ij 
and Sq^q' are the Kronecker symbols. Properties (i) and 
(ii) follows from the Perron-Frobenius theorem, see, for 
example, Ref. [11]. Properties (iii) and (iv) are proved 
in Appendix |^ One can show that for an uncorrelated 
network there is a single largest eigenvalue and an {Nd — 
l)-degenerate zero eigenvalue: 

Ai = Z2/Z1, Xi>2 = 0. (16) 

The percolation threshold corresponds to a critical prob- 
ability Pc above which a nontrivial solution of Eq. (jl2p , 
Vq > 0, appears. Taking into account only the linear 
terms, we get the following condition: 

yq^P^Bqq>yq.. (17) 
9' 

We represent yq as a linear combination of the mutually 
orthogonal eigenvectors which we call modes: 

i=l 



The amplitudes are unknown functions of p. It is 
obvious that Oj = at p < Pc- Substituting Eq. (fTS]) 
into Eq. (fTT]) and using the orthogonality Eq. (|14p of the 
eigenvectors, we obtain an equation a.i = pXiUi. One 
can see that a nontrivial solution ai ^ appears when 
pXi = 1. So the mode associated with Ai is critical. This 
gives the following criterion for the percolation threshold 
found in Refs. [S, .16]: 

PcXi = 1. (19) 

Thus the generalization of the MoUoy-Reed criterion to 
undamaged correlated networks, i.e., at p — 1, is the 
following condition: if the largest eigenvalue Ai o/ the 
branching matrix, Bqqt , is larger than 1, then the corre- 
lated network has a giant connected component. In un- 
correlated networks the criterion (fT9|) is reduced to the 
well-known one: Z2I z\ — 1. 

In Appendix we prove that at a given degree distri- 
bution P((7), the largest eigenvalue A^'^'''' of an assortative 
network is larger then z^j z\ while in a disassortative net- 
work it is smaller: 

Af''<-<Aj"^\ (20) 
^1 

There are the following lower and upper boundaries for 
Ai (see Appendix E]) : 

- VB(q)q(q-l)P(q) < Ai <maxB(g). (21) 

Suppose that a degree distribution is such that Z2I z\ — 1. 
Then inequality (pO)l leads to the following statements, 
(i) An uncorrelated network with this degree distribution 
is at the birth point of a giant connected component, (ii) 
An assortative network with this degree distribution has 
a giant connected component, (iii) A disassortative net- 
work with this degree distribution has no giant connected 
component. According to criterion pop and inequality 
(1^01) , percolation thresholds in randomly damaged assor- 
tative and disassortative networks with the same degree 
distribution satisfy the inequality: 

pi-^^) < zi/z2 < (22) 

Therefore, assortative mixing enhances resilience of a cor- 
related network against random damage while a disassor- 
tative mixing decreases it [9[. 

If the largest eigenvalue Ai diverges in the infinite net- 
work, then the percolation threshold p^ tends to zero. 
In this case the giant connected component of a cor- 
related network cannot be eliminated by a random re- 
moval of vertices p^ . In an assortative network, in ac- 
cordance with inequality (PO)) . this takes place when the 
second moment Z2 diverges. This criterion (z2 — > 00) of 
the robustness of an uncorrelated network was found in 
[13, [13] ■ Remarkably, Ai can diverge even if an assor- 
tative network has a finite Z2 but degree-degree correla- 
tions are sufficiently strong. An example of an assorta- 
tive network of this kind is given in Sec.|Vl On the other 
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hand, Vaazquez and Moreno ^16*1 found a network with 
strong disassortative mixing, which has a finite percola- 
tion threshold and is fragile against a random damage 
even if the second moment zi diverges. 



IV. CRITICAL BEHAVIOR OF "WEAKLY" 
CORRELATED NETWORKS 

In this section we derive necessary and sufficient con- 
ditions for degree-degree correlations to be irrelevant for 
critical singularities at the percolation transition. These 
conditions are: 

(I) The largest eigenvalue Ai of the branching matrix 
must be finite if (g^) is finite, or Ai ^ oo if (g^) 
oo. 

(II) The second largest eigenvalue A2 must be finite. 

(Ill) A sequence of entries ^q^'' of the maximum eigen- 
vector of this matrix must converge to a finite non- 
zero value at 5 — > 00, see Eq. (|30l) . 

Let us solve Eq. near the percolation transition. 
We use the fact that the eigenvectors $g*'' of the branch- 
ing matrix B form the complete orthogonal basis set, see 
Eqs. dUl) and nil). Substituting Eq. HH) into Eq. p2)) . 
we get a set of nonlinear equations for unknown ampli- 
tudes ai: 

{pX,-l)a,=pX,S,{a). (23) 
Here Si{a.) is a function of the amplitudes a — (ai, a2, ...): 

Zl 



q=3 n=2 



1 



{-y,rqP{q)^^\ (24) 



Let us find a; in the leading order in r 1 under con- 
ditions I and II. Substituting Eq. (HI]) into Eqs. ((221) and 
(|24p and taking into account only the quadratic terms in 
tti , we get a set of approximate equations for the ampli- 
tudes a,-: 



[pXi - l)ai = pXi Mirnnaman, (25) 



where 



M^^n = ^ E '^(9 - 1)(9 - 2)P(g)$«$(")$ 



(m) 
1 ■ 



(26) 



Let us first consider a degree-degree correlated network 
with finite coefficients Mimn in the infinite size limit N 
oo. In the leading order in r, equation ()25|1 has a solution 



One can see that the amplitude ai of the critical mode 
Ai has the standard mean-field dependence. It is much 
larger than the amplitudes ai of the modes i > 2 which 
we call transverse modes. Therefore near the percolation 
threshold the order parameter = 1 — Xq is mainly 

determined by the critical mode ^'"q^ associated with the 
largest eigenvalue. The transverse modes give a smaller 
contribution: 



ai$W +0(a?). 



(29) 



Above we assumed that the coefficient Mm is finite. 
This assumption takes place if the sequence of entries 



$5^^ converges at q 



qcut 



cxD, e.i., lim„ 



< oo. 



If this limiting value is non-zero, 
< lim < 

q — *oo ^ 



(30) 



then the third moment {q^) of the degree distribution 
must be finite, (g^) < co. Equation ([30|) is the condi- 
tion III formulated above. Thus we conclude that under 
conditions I-III the percolation in a correlated scale-free 
networks with 7 > 4 has the standard mean-field critical 
behavior, Eq. (|77|) . with the critical exponent (3 = 1. In 

Sec. IVIII we will show that the case limq^oo '^'q^^ — 0, 
e.i., the condition III is broken down, corresponds to a 
strongly correlated network. 

Consider the case 3 < 7 < 4 under conditions I-III. In 
this case it is necessary to take into account all orders in 
ai in Eq. (|23p . The reason is that the coefficients of the 
expansion over ai diverge as A'^ — > 00. It is this diver- 
gence that leads to the non-standard critical behavior of 
percolation on an uncorrelated complex network [Tsj . 
We assume that the inequality ai ^ for alH > 2 is also 
valid at 3 < 7 < 4. Below we will confirm this assump- 
tion. We expand the function S'i(a) over the amplitudes 
ai of the transverse modes and save only linear terms in 



where 



S'i(a) ~ Si(ai) + ^ Sij{ai)aj 

3>2 



Si(ai) = Si{ai,0,0, . . .), 
_ dSiia) 

Stj[ai) = la=(ai,0,0,..)- 



(31) 

(32) 
(33) 



The function Si(x) is given by the following series: 

sdx)= K,{n)(~xr, (34) 



ai 



Mm 
XialM,i 



(27) 



(Ai - AO 



i cx 0{t^) for i > 2. (28) 



where 



Kdn) = - J2 ('^ ^)gP(g)($W)"cI>W. (35) 



q=n+l 



6 



Similarly, the function Sij(x) is 



J2 K,,{n){-x)\ 

n=2 



(36) 



At 3 < 7 < 4 the functions Si (x) and Sij (x) are singular 
functions of x because the coefficients Ki (n) and Kij (n) 
diverge as cicut(N) — » oo. Using the method of summa- 
tion used in [ij, [3^ we get asymptotic results: 



^(x) ~ biX^ , 



(38) 



Here the numerical coefficients and 6^, where i,j = 
1,2, ... , depend on a specific model of a correlated net- 
work. Substituting Eqs. ^ and §^ into Eq. ^ we 
obtain a set of nonlinear equations 



{pXi - l)a.i pA^foia? ^ + 



,7-3 



J>2 



for i — 1, 2, ... . They have an asymptotic solution 



1/(7-3) 



ai 



"(aT^X) 



oc 0{t 



(7-2)7(7-3) 



(39) 

(40) 
(41) 



In the case 3 < 7 < 4 the amplitude ai of the critical 
mode ^[^^ also is mush larger than the amplitudes of the 
transverse modes, i.e., ai ^ ai for all i > 2. So, in the 
leading order, we get yq « ai^[^\ The critical exponent 
P is equal to 1/(7 — 3) as for an uncorrelated network 
with the identical degree distribution. 

Now consider the case 2 < 7 < 3 under conditions (I)- 
(III). According to Eq. ((20|) . in an assortative network 
with a divergent second moment Z2 the largest eigen- 
value Ai diverges in the limit N ^ 00. The percolation 
threshold is Pc = 0, and a giant connected component 
is present at any p > 0. One can show that at p <C 1 
the order parameter i/q and the size of giant connected 
component S are 



S oc p^+^ 



(42) 
(43) 



where (3 = 1/(3 — 7). This is the same singularity as was 
found in uncorrelated networks with 2 < 7 < 3 [H, [lB| ■ 
This result is also valid for disassortative networks with 
Ai = 00. 

The main conclusion of this section is that under con- 
ditions (I)-(III), a correlated complex network has the 



same critical behavior as an uncorrelated random net- 
work with the identical degree distribution. Our results 
are summed up in Table HI 

Degree-degree correlations may be strong and violate 
one of these conditions. In Sec. IVII and Sec. IVIll we will 
give examples of such networks and new critical singu- 
larities. 



V. EXACTLY SOLVABLE MODEL OF A 
CORRELATED NETWORK 

In this section we introduce a simple model of a cor- 
related network which allows analytical treatment. This, 
actually, toy model permits us to check general results 
of the percolation theory derived above and to study an 
effect of strong assortative mixing. This model is also 
interesting due to its clear modular structure. 

Let us consider a correlated network with a given de- 
gree distribution P{q) and a degree-degree distribution 
P{q,q') which is factorized a.t q ^ q' as follows; 



P{q,q')^Sq,q,[f{q) + l]p'iq)P'{q) 
+ il-Sq,q')p{q)Piq)p{q')P{q'). 



(44) 



Compare this with an uncorrelated network where 
P{q,q') — qP{q)q' P{q')/ z1 which corresponds to p{q) = 
qj z\. The function /(q) is non-negative for assortative 
mixing. In this case vertices at the ends of an edge have 
the same degree with a higher probability rather than 
different degrees. In the limit /(g) 00 we have a set 
of disconnected random g-regular modules. For disas- 
sortative mixing, we choose —1 < /(q) < 0, then the 
probability to have different degrees at the ends of an 
edge is higher than to have the same degree. 

For a given degree distribution P(g) and a function 
!{q), the substitution of Eq. dM]) into Eq. (P) gives the 
following equation for the function p(q): 



p{q) 



zipl + [l + 4g/(g)P(g)/(zip)]i/2 



(45) 



where the parameter p must be found self-consistently 
from the equation: 



= Y^p{q)P{q). 



(46) 



For model (|44)) . eigenvalues A and eigenvectors ^')^^ of 
the branching matrix, Eq. Ill, are determined by the 
following exact equations: 



{q-\)p\q)P[q) 
q A-A(g)^ 



^(A) ^ zxp{q)n{\) 

" g[A-A(g)]- 



(47) 
(48) 
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FIG. 2: Dependence of the largest, Ai, and the second largest, 
A2, eigenvalues of the branching matrix B on the parameter / 
in the uniform case f{q) — f for model (|44p with a scale- free 
degree distribution P{q) = kg""*, {q) = 2.3, {g(g- 1)) = 11.3. 
Other eigenvalues are non-negative and lie below A2. The 
eigenvalues are normalized in respect to Ao = z-ijzx. 

Here we have introduced a function 

A(g) = ii^zi/(q)p2(q)p(^). (49) 

n(A) is a normalization constant determined by Eq. ()14p . 
Choosing a function /(g) and degree distribution P{q), 
we obtain a correlated network with a specific spectrum 
of the eigenvalues of the branching matrix. 

For model (|44|) the average branching parameter B{q) 
defined by Eq. Q is equal to 

Big) = ^iM ^(5' _ l)p(g')P(g') + A(g). (50) 

Assuming that P{q)f{q) ^ 1 at large q, we obtain the 
asymptotic behavior: p{q) « q/{pzi). Therefore, B{q) — 
const -I- A(g) at g > 1. If A(g) cx and -1 < a < 0, 
B{q) increases with increasing q. 

Model (|44|) has a modular structure. Modules are 
formed by vertices with equal degrees q. The number of 
modules is equal to the number of different degrees, Nd- 
In the limit /(g) — > 00 for all g, this network consists of 
disconnected modules. Each module is a g-regular ran- 
dom graph. In this limit, the Pearson coefficient r = 1, 
which corresponds to a perfectly assortative network Q . 
If /(g) is large but finite we get a weakly connected mod- 
ules. 

Equations (|T7|) and show that the largest eigen- 
value Ai is larger than maxA(g). All other eigenvalues 
lie in the range: minA(g) < A < maxA(g). For ex- 
ample, if the function A(g) is a monotonously increasing 
function, then there is only one eigenvalue in every inter- 
val (A(gi), A(gi-|_i)). If maxA(g) < 00, then A2 is finite 



FIG. 3: Relative size S of a giant connected component versus 
the occupation probability p for correlated network (|44|) hav- 
ing two groups of vertices with degrees qi — 4 and 92 ~ 11, 
P{qi) ~ 0.6 and P(q2) = 0.4. The curves were obtained 
by numerical solution of Eq. (|10|l in the following cases: (i) 
/(91) = fil'i) = 0, i.e., an uncorrelated network (solid line); 
(ii) an assortative network with /(gi) = /(g2) = 10 (dotted 
line); (iii) an assortative network with /(gi) = fil'i) ~ 100 
(dashed line). Arrows indicate the positions of I/A2 for (ii) 
and (iii). 




FIG. 4: Amplitudes ai and 02 of critical and transverse 
modes, respectively, versus p — pc in an assortative network 
with degree-degree distribution (|44p and the same parame- 
ters as for the dotted curve in Fig. [3] f{qi) ~ f{q2) ~ 0. 
The percolation threshold is Pc = 1/Ai . The arrow shows the 
point p = I/A2. 



because A2 < maxA(g). Ai is finite at (g^) < 00 and 
diverges at (g^) — > 00. Equation (|48|) shows that the 
maximal eigenvector satisfies Eq. pO|) . Thus conditions 
I-III in Sec. IIVI are fulfilled if maxA(g) is finite. In this 
case model (|44p belongs to the same universality class as 
an uncorrelated random complex network with the iden- 
tical P{q), e.i., degree-degree correlations does not affect 
the critical behavior. 
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In Fig. [2] we present the results of numerical calcu- 
lation of the largest and the second largest eigenvalues 
of the branching matrix for a scale-free degree distri- 
bution, P{q) — Aq~'^, and uniform assortative mixing 
fio) — f > ^- One can see that Ai and X2 are finite even 
correlations are strong, f{q) = f^l- 

In order to check our predictions, Eqs. (|27p and (pS)) . 
for the amplitudes of critical and transverse modes at the 
percolation threshold, we solved numerically Eq. (flO)) for 
model Eq. We considered a simple network con- 

sisting of two groups of vertices with degrees qi = 4 and 
(72 = 11 and the following degree distribution: P{qi) = 
0.6, P{q2) = 0.4. In this situation the branching matrix 
has only two modes — a critical mode and a transverse 
one. The resulting dependence of the size S of a giant 
connected component on p is shown in Fig. [3] A giant 
connected component emerges at pc = l/^i. One can 
also see a non- monotonous behavior of dS{p)/dp which 
takes place when the occupation probability p crosses 
I/A2. Figure m shows that the amplitudes of the critical 
and transverse modes near the percolation point behave 
in agreement with Eqs. (|77|) and (|^ . i.e., ai oc {p — Pc) 
and a2 oc (p — Pc)^. 

VI. CRITICAL BEHAVIOR OF STRONGLY 
ASSORTATIVE NETWORKS 



0.5 





1 




2 /^^^ 









FIG. 5: Phase diagram in the plane (7,0) of model (|44|l . 
Here 7 is the degree distribution exponent, P{q) (x g""*", and 
a is the exponent of the function A{q) oc (?~". In region I 
the critical behavior of the model is the same as that for an 
uncorrelated network with the identical degree distribution. 
In region II the sequence of eigenvalues is unbounded, and this 
model demonstrates a new critical exponents, see TableU] The 
boundary between regions I and II is defined by the equation 
a = min(0, 7 — 8). 



Let US consider the case when the function A{q) is un- 
bounded, i.e., maxA(g) — 00. Then the sequence of 
eigenvalues Ai also is unbounded. In particular, both 
Ai and A2 are infinite. The condition II is not fulfilled. 
For example, we can choose the function A(q) = Aoq^". 
It is unbound at — 1 < a < and bounded at a > 0. In 
the latter case, model (|33]) has the same critical behav- 
ior as an uncorrelated network with the identical degree 
distribution. In Fig. [5] this region of parameters a and 
the degree distribution exponent 7 is defined as region I. 

In order to study percolation in model (jH]) for — 1 < 
a < 0, we rewrite Eq. (jlOp in the following form: 



pA{q) r 



q-1 



1- l-?^^ 



zip{q)v 



(51) 



where the parameter v plays the role of an effective field: 



■pY^p{q)P{q) l-[l-y. 



E 



qvq 



zip{q)f{q)' 



/(?) 



(52) 



At p <JC 1 and qcut 2> 1 these equations have an approxi- 
mate solution 

pA{q) zip{q)v 
Vq = T H , at g > gp. 



q-l q 

zip{q)v 
q[l-pA{q)V 



at g qp, 



(53) 



where a characteristic degree qp is determined from the 
equation pA(qp) = 1 which gives 



= (Aop)!/" » 1. 



(54) 



Within this solution we have {q — l)yq S> 1 at q 3> qp, 
though <C 1. This enables us to use the approxima- 
tion: 



1 - 



q-l 



1. 



(55) 



For the degree distribution P{q) = Aq with 7 > 3, a 
self-consistent solution of Eq. (|52p is 



Ap 



ziPil - '2)qp 



,7-2 



oc p 



(2-7+q)/q 



(56) 



The leading contribution to v is given by a sum over 
degrees q > qp. Substituting solution ([53|) into Eq. (jlip . 
we obtain the size of a giant connected component: 



(57) 



Interestingly, the main contribution to 5* comes from de- 
grees q < qp. 

For a slowly increasing function A(g) ~ Aq In q (which 
corresponds to a 0), we have qp = exp(l/(pAo). Then, 
at an arbitrary 7 > 3, we obtain from Eq. (|56p the fol- 
lowing results: 



wocpe-(^'2)/(pAo)^ 
5ocpwocp2e-(T-2)/(pAo)_ 



(58) 
(59) 
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On the other hand, for 2 < 7 < 3 we note that if qcut 
00, then the series Eq. (|52p is a singular function of v. 
Then, Eq. takes the form: 

V « apv^'^ + . ^_2 ^ (60) 

where a = 0{1) is a model dependent parameter. Here 
the first singular term is given by the sum over q < qp 
under the condition that qpV <^ 1. The second term is 
given by degrees q > qp. At7>3 + a we arrive at the 
solution Eq. ([56]) while at2<7<3 + awe obtain 

V (X (61) 

with (3 = 1/(3 — 7). The size of a giant connected com- 
ponent is S" oc oc p^~^^ . This is exactly the behavior 
that was found for an uncorrelated scale-free network for 
2 < 7 < 3 [l3, [3 ■ It means that degree-degree correla- 
tions are irrelevant at2<7<3-|-Q!. 

Our results are summed up in Fig. [5] and Table IH In 
region I, degree-degree correlations are irrelevant. They 
become relevant and lead to new critical exponents in the 
region of parameters 7 > 3-1- a and —1 < a < (region II 
in Fig. [5]) where model ([ii]) has an unbounded sequence 
of eigenvalues of the branching matrix. 



1 n 




2 3 4 



Y 



FIG. 6: Phase diagram in the plane (7, a) of model (|62|l . Here 
7 is the degree distribution exponent, P{q) oc q and a > 
is the exponent of the function g{q) oc g~". In region I this 
model shows the standard mean-field critical singularity with 
the critical exponent [3 = 1. In regions II and III this model 
demonstrates a new critical behavior with the critical expo- 
nent depending on a, see Table IT] The boundaries between 
regions I and II, and II and III, are defined by the equations 
Of = (4 — 7)/3 and q = (3 — 7)/2, respectively. 



VII. CRITICAL BEHAVIOR OF STRONGLY 
DISASSORTATIVE NETWORKS 

Vazquez and Moreno introduced networks with the 
following degree-degree distribution: 

P{q',q) = -(1 - gq')q'P{q'){l ~ 
^1 



-{l-g,)qP{q){l-5,.j)5q,.j 



1 

zfG 



-^g,-q'P{q')g,qP{q){l 



V,i)(l-'5,,l), (62) 



where G = ^gyi gqqP{q)/ zi. The fraction of vertices 
of degree 1 is fixed by the following condition: P{q — 
1) = X]<j>i(l ~ 9q)<}P{Q)/ zi- In this model, a vertex 
of degree q > 1 is connected to a vertex of a degree 
q' — 1 with probability 1 — g^, and this vertex (q > 1) is 
connected to a vertex of degree q' > 1 with probability 
9q9'q I G. From Eq. ^ we obtain that the average nearest- 
neighbor's degree of the vertices of degree q is B{q) oc 
gq. These networks are disassortative for any monotonic 
decreasing function gq. 

One can show that the branching matrix has the 
largest eigenvalue 



ZlLz 



(63) 



All other eigenvalues are zero, Ai = for i > 2. Remark- 
ably, this spectrum is similar to that of an uncorrelated 



network, see Eq. (|16p . The maximal eigenvector is 



$(1) = 



(64) 



where ni is for normalization. Using this result, we find 
the local clustering coefficient C{q), Eq. (|B3|) . of model 
(EH): 



Ciq) 



A? 



NzfG 



(65) 



Let us study the case of a scale-free degree distribution 
P{q) oc and choose 



9q 



(66) 



with a > 0. Zero exponent a = corresponds to an 
uncorrelated network. If a > 0, then — > at g ^ 00. 
It means that condition III in Sec. IIVI is not fulfilled. 

The dependence of the percolation threshold Pc on the 
parameters 7 and a was found in ^6]. Let us study the 
critical behavior of model ((62)) . From Eq. ([23]) we find 
that the amplitudes = for i > 2. Therefore, the 
order parameter is completely determined by the critical 
mode: = ai^'^\ The amplitude ai is determined by 
the following equation: 



(pAi 



Dai 



pAiSi(oi), 



(67) 



where the function si{x) is given by series Eq. p4p . First 
we study the region of parameters a and 7, where 



a > max(0, ■ 



■1. 



7 > 2 



(68) 
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(region I in Fig. [6]). Based on Eqs. (|26p and (|63|) one 
can see that in this region, the coefficient Mm and Ai 
are finite. Therefore, the percolation threshold is non- 
zero, Pc = 1/Ai ^ 0, e.i., this network is fragile against 
a random damage. The Eq. (j67[) has a solution oi oc 
with standard mean- field exponent P — 1. Remarkably, if 
a > 2/3, then this standard critical behavior takes place 
at all 7 > 2. Region I includes a subregion with 2 < 7 < 3 
where the second moment (g^) of the degree distribution 
diverges. This is in contrast to an uncorrelated network 
with 2 < 7 < 3, which has Pc ~ and exponent (3 — 
1/(3-7). 

Region II in Fig. [H] is defined by the inequalities 



7 4 - 

— - < a < — 



■ 7 



2 < 7 < 4. 



(69) 



In region II, Ai is finite while Aim — 00. This points 
out that the function si{x) is a singular function. The 
asymptotic behavior of this function is the following: 
si{x) oc x^^" where = (4 - 3a - 7)/(l - a) > 0. 
Equation (j67[) gives oi oc with the exponent 



f3 



1 — a 



j-3 + 2a 

Region III in Fig. [6] is defined by the inequalities 
3-7 



< a < 



2 < 7 < 3. 



(70) 



(71) 



In this region both Ai and Mm diverge. Therefore, Pc — 
and the network is robust against random damage. To 
find the order parameter we directly solve Eq. pO|) . The 
exact solution is i/q = ygq where y satisfy 



y 



P 
ziG 



?>i 



(72) 



The series in y on the right-hand side has an asymptotic 
behavior oc y^~^^^ ^ y <C 1, where 



1 



3 - 2a - 7 



(73) 



This leads to a solution y (x p^^ . The resulting size of a 
giant connected component is S* oc py. 

Thus, strong disassortative correlations can dramati- 
cally change both the percolation threshold and the crit- 
ical behavior. In a broad region of parameters 7 and 
a, see Fig. Inland Tabled the critical properties of model 
(I62p with disassortative mixing differ from those of a cor- 
responding uncorrelated network. 



VIII. CONCLUSIONS 

We have studied critical phenomena at the percolation 
transition in equilibrium complex networks with degree- 
degree correlations. Our consideration is based on the 



assumptions that a network is locally tree-like, i.e., clus- 
tering is negligibly small, and that only degree-degree 
correlations between the nearest neighboring vertices are 
present. The origin of degree-degree correlations is not 
relevant for our approach. These correlations can be in- 
trinsic (i.e., they implicitly follow from other features of 
a network) or directly defined in a given network model. 

We have demonstrated that both assortative and disas- 
sortative mixing affect not only the percolation threshold 
but can also change critical behavior at this percolation 
point. We have found necessary and sufficient conditions 
for a correlated network to have the same critical be- 
havior as an uncorrelated network with the identical de- 
gree distribution. These conditions result from the fact 
that critical singularities at the percolation point are de- 
termined not only by a degree distribution, as for un- 
correlated networks, but also spectral properties of the 
branching matrix. The resulting critical behavior of a 
correlated network belongs to the same universality class 
as percolation on an uncorrelated network if the follow- 
ing conditions are fulfilled: (I) the largest eigenvalue Ai of 
the branching matrix is finite if (q^) is finite, or Ai — > 00 
if {q^) 00; (II) the second largest eigenvalues of the 
branching matrix is finite; (III) the sequence of entries 
of the eigenvector associated with the largest eigenvalue 
converges to a non-zero value. In this situation, one can 
say that degree-degree correlations (assortative or dis- 
assortative) are irrelevant for critical phenomena though 
they change the value of a percolation threshold. The 
critical exponents are completely determined by asymp- 
totic behavior of degree distribution at large degrees. 

Degree-degree correlations become relevant if at least 
one of these conditions is not fulfilled. We have pro- 
posed two simple models of correlated networks with 
strong assortative and disassortative mixing, where cor- 
relations dramatically change the critical behavior. As 
a result, the critical exponent becomes model dependent 
and hence non-universal. The advantage of these mod- 
els is that they allow us, first, to change gradually from 
weak to strong mixing, and, second, to obtain the ex- 
act analytical solution of the percolation problem. In 
the assortative networks, proposed in this work, strong 
assortative mixing leads to an unbounded sequence of 
eigenvalues of the branching matrix, e.i., condition II is 
not fulfilled, which results in new critical singularities. 
Remarkably, in these networks the percolation threshold 
can be zero despite a finite second moment (g^) of a de- 
gree distribution, in contrast to uncorrelated networks. 
(In uncorrelated networks, Pc = only if (g^) — > 00.) 
We have also found an unusual critical behavior in a net- 
work with strong disassortative mixing where condition 
III is not fulfilled. Here the situation is just opposite 
to the strongly assortative network. Namely, in contrast 
to uncorrelated networks, the considered disassortative 
network can be fragile against random damage even if 
(q^) 00. Moreover, we have shown that in a wide 
range of model parameters, strong disassortative mixing 
results in new critical singularities. 
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Many real-world networks and network models have 
non-zero clustering and long-range degree correlations. 
Intrinsic degree correlations are present, for example, in 
the static model of scale-free networks with degree expo- 
nent 2 < 7 < 3 [Slj and in some other models without 
self-loops and multiple connections. Unfortunately, these 
networks are clustered, and it is unknown whether their 
degree-degree correlations are short- or long-ranged. A 
generalization of the theory to networks with long-range 
degree correlations and non-zero clustering is a challeng- 
ing problem [Toj . 

We believe that our general conclusions concerning a 
strong influence of degree-degree correlations on criti- 
cal singularities are true for various models of statistical 
physics defined on the top of correlated networks. 
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APPENDIX A: SPECTRAL PROPERTIES OF 
THE BRANCHING MATRIX 

Here we describe general properties of the eigenval- 
ues Xi and eigenvectors <i>^*'' of the branching matrix ([4]). 
According to the Perron-Frobenius theorem, a real posi- 
tive matrix B has a positive real eigenvalue Ai such that 
Ai > |Ai| for any eigenvalue Ai of B, see, e.g., Ref. [36| . 
Furthermore, there is an eigenvector with positive entries 
^[^^ > corresponding to Ai. Ai is the largest eigenvalue. 
The eigenvector associated with Ai is called a maximal 
eigenvector. If this matrix is positive then the largest 
eigenvalue is simple, i.e., non-generate. 

One can show that a symmetric matrix 



Bq,q' - 



(Al) 



has the same eigenvalues A^ as the non-symmetric matrix 



B, i.e., Bqq'^q' = ^i^q \ whcrc eigenvectors $q are 



(i) 



related to the eigenvectors $q as follows: 

$W = [g(g-l)P(q)/(q)]l/2$W. 



(A2) 



This relationship between the non-symmetric matrix B 
and the symmetric matrix B allows us to obtain the 
following spectral properties. First, since all eigenval- 
ues of a symmetric real matrix are real, all eigenvalues 
Ai of B are also real. They can be ordered as follows: 
Al > A2 > ... > Ajv^. Second, the eigenvectors <I>'^*) form 



a complete orthonormal basis set: 
q 

^q ^q' 



(A3) 
(A4) 



Using relation (|A2p we obtain Eqs. and (HSl). Be- 
cause the maximal eigenvalue Ai of the positive matrix B 
is non-degenerate, this eigenvalue is strictly larger than 
the second largest eigenvalue A2. In other words, there is 
a gap between Ai and the second largest eigenvalue A2, 
Al — A2 > 0. A value of this gap is determined by a spe- 
cific form of a joint degree-degree distribution P(g', q). 

The largest eigenvalue Ai of a non-negative branching 
matrix Bq^i has the following upper boundary [36| : 



Al < max_B((7). 
q 

Using the equation 

Eqq' ^q'{q'-l)P{q\q){q~'l)^q 



(A5) 



(A6) 



which directly follows from Eq. (|13p . we can find a lower 
boundary. Substituting 'J, = 1, we get 



1 



Z2 



B{q)q{q - l)P{q) < Ai 



(A7) 



Unfortunately, no results are known for the second largest 
eigenvalue A2. 

For an uncorrelated network, P{q' , q) — 
q' P{q')qP{q)/ {q)'^ . In this case the matrix B has 
the entries 

B,,, = {q'^l)q'P{q')/{q)^B^^;^,. (A8) 
The largest eigenvalue of this matrix is 



Al = Z2/Z1. 



(A9) 



The normalized maximal eigenvector has equal entries 



/2 



(AlO) 



for all q. There is also an {Nd — l)-degenerate zero eigen- 
value, A2 = ... = \ncl = 0- In this case the gap, Ai— A2, in 
the spectrum is equal to 22/^1 ■ Any vector <I> orthogonal 
to in the sense of Eq. ^ is an eigenvector associated 
with this zero eigenvalue. 

Let us consider a network in which a joint degree- 
degree distribution P{q',q) slightly deviates from 
q'P{q')qP{q)/{qf 



I.e. 



max ■ 

q 



\B{q)~B\ 



B 



< 1, 



(All) 



where the branching coefficient B{q) and the mean 
branching coefficient B are defined by Eqs. (O and ([5]). 
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In order to find the largest eigenvalue, we write B = 
B^°'> + 6B, where the matrix B^°^ is given by Eq. (|X8|) . 
and 613 = B — B'^'^^ is a perturbation. In the first order 
of the perturbation theory we get: 

$W « AB{q)/B, (A12) 
X,^-Y^B{q)q{q-l)P{q). (A13) 



Here A is a normalization constant. Equation (|A13[) 
agrees with the lower boundary in Eq. (|A7|) . Let us ana- 
lyze these results. First, the entries <^'^(p of the maximal 
eigenvector are positive and lie in a bounded range. Sec- 
ond, we can rewrite the right-hand side of Eq. (jA13p as 



X.^B^ = 



Y,{q' ~l)P{q'^q){q^l)-B' 



= B + ra^/B. 



(A14) 



We see that degree-degree correlations give a contribu- 
tion which is proportional to the Pearson coefficient r, 
Eq. ([9]). Therefore, assortative mixing increases Ai while 
disassortative mixing decreases it. As a result we get 
inequality (I20p . With increasing degree-degree correla- 
tions the magnitude of the second largest eigenvalue A2 
increases gradually from zero. Therefore, in a weakly cor- 
related network the matrix B has a finite gap Ai — A2 > 
which depends on a specific degree-degree distribution 



APPENDIX B: RELATIONSHIP BETWEEN THE 
MATRIX B AND NETWORK PARAMETERS 

Many structural parameters of correlated networks can 
be related with the spectrum of the branching matrix. 
The local clustering coefficient C{q) of a vertex with de- 
gree q is defined as follows: C{q) = t{q)/[q{q — l)/2], 
where t{q) is the number of triangles (loops of length 3) 
attached to this vertex, and q{q — l)/2 is the maximum 
possible number of such triangles. The clustering coef- 
ficient C and the mean clustering C are given by the 
following equations: 



C^-Y.{q-l)P{q)C{ql 



where B is the branching matrix ([4]) . Using this equation 
and Eq. (|15p , we find relationships between the clustering 

coefficients and the eigenvalues A^ and eigenvectors $g 
of the branching matrix: 



N{q)^^ 



1 



i—1 



(B3) 



(B4) 



C 



1 



^A?^P(g)($«)^ (B5) 



For an uncorrelatcd network, C^"") = z^/{N{qf) 
In a network with assortative mixing the largest eigen- 
value Ai is larger than 22/2:1 according to Eq. (|20p . 
Therefore the clustering coefficient C'^"") is larger than 
(j(un) ^ Similarly, in networks with disassortative mix- 
ing the clustering coefficient C^'**^-' is smaller then C'-""-'. 
Thus we obtain an inequality: 



< 



N{qy 



< 



(B6) 



Using Eq. IjlSp. we obtain a useful relation between 
the average branching coefficient B{q), Eq. ([6]), and the 
eigenvalues A^ and eigenfunctions $, 



(B7) 



where Vi = ^qW[q)^q\ This relation shows that the 
degree dependence of B(q) is determined by the eigen- 
vectors mainly by the maximal eigenvector $5^^ 
Note that if Ai 00, then B{q ~> 00) 00 but not vice 
versa. 

The average number, Ni(q'\q), of vertices with degree 
q' at a distance I from a vertex of degree q can also be 
related with the branching matrix B: 



Nl{q'\q) 



Q'-l 



(B8) 



The mean intervertex distance £ of the correlated network 
can be found from the equation: 

1 

q' n—1 



C = Y,P{q)C{q). 



(Bl) 



According to Ref. [40| , we can rewrite C [q] in the follow- 
ing form: 



C{q) 



Nq{q-l)Piq)- 



(B2) 



where 



A n=l 



(B9) 



E 9(9 -1)^(9)* 



(BIO) 
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In a network with a finite gap between the largest the 
second largest eigenvalues, A2, the leading contribution 
into the right-hand side is given by the terms with Ai. 
As a result we obtain the following estimate of 1 for a 
correlated network: 



In TV 

hTvT' 



(Bll) 



Note that an uncorrelated network has a mean interver- 
tex distance £ w InN/ ln(z2/zi). Using (j20|) we find that 
for a given degree distribution P{q) a network with as- 

sortative mixing has a smaller intervertex distance 1 
than an uncorrelated network, while a network with dis- 

assortative mixing has a larger £^ ' : 



jias) ^ 



InN 



ln(z2/zi) 



<£ 



(di 



(B12) 



APPENDIX C: ENTROPY OF A CORRELATED 
NETWORK WITH A MODULAR STRUCTURE 

Here we describe a statistical ensemble of assortative 
networks with degree-degree correlations, Eq. (|44)) . The 
probability Pa(ay) that an edge between vertices i and 
j of degrees qi and qj is present (oy — 1) or absent 
(fly = 0) is 



given. Then the probability of realization of a graph with 
a given adjacency matrix a^- is the product of probabili- 
ties Pa(oy) over all pairs of vertices: 



N-l N N 

^(K}) = ^n n PMX{5{^a,, - q^. 

1=1 i=i+i 1=1 i 

(C2) 

The delta-functions fix degrees of the vertices. Z is a 
normalization factor: 



Z^exp[^E^(^)l4^(^) 



(C3) 



where g - {q)[l + p'' + Y.q -P' (<?)/(«)]■ In fact, Z is 
the partition function of this network ensemble (compare 
with a partition function for other network ensembles in 
Ref. [3§|). The function p{q) is given by Eq. (gS]). At a 
given (g), parameters p and /(g) are coupled. The con- 
dition of the maximum of the entropy In Z with respect 
to the parameter p leads to Eq. 

The average of a physical quantity A[{aij}) over the 
network ensemble is given by 



N-1 N 



A{{a.,,})V{{a.,,}) n n (C4) 
2—1 j—i-yi 



Zip 
N 
2 



(l + /(%)'5,.,,,)^(ay-l) 



Zip 
N 



(CI) 



Here are the entries of the adjacency matrix, p and 
f{q) are parameters. For an uncorrelated network (the 
configuration model), we have p = 1 and /((?) = 0. A 
positive f{q) results in assortative mixing. Let a se- 
quence of degrees [and so a degree distribution P{q)] be 



In particular, for a graph with a given adjacency matrix 
Qij, the degree-degree distribution P{q, q') is given by an 
average over edges. 



Piq,q') ^ J^J2^1^'1'''i^1^'1'■ 



iC5) 



Averaging this function over the network ensemble, we 
get Eq. 
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